Arbitrary amplitude inertial Alfven waves in 
homogeneous magnetized electron-positron-ion 

plasmas 

S. Mahmood and H. Saleem 
Department of Physics, COMSATS Institute of Information Technology, 

H-8 Islamabad Pakistan 
and 

PINSTECH (Physics Research Division) P.O. Nilore 
Islamabad Pakistan 

February 2, 2008 

Abstract 

Nonlinear set of equations for inertial or slow shear Alfven wave (SSAW) 
in ideal electron-positron-ion (e-p-i) plasmas are presented. The analytical 
solution for arbitrary amplitude SSAW in such multi-component plasmas 
is obtained using Sagdeev potential approach. The numerical solutions for 
several different cases have also been presented for illustrative purpose. 
It is found that the electron density dips of SSAW are formed in the su- 
per Alfvenic region. The amplitude and the width of the nonlinear shear 
Alfven wave reduces with the increase in the concentration of positrons in 
electron-ion (e-i) plasmas. The width of the soliton also depends upon the 
direction of propagation of the perturbation in both e-i and e-p-i plasmas. 
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I. INTRODUCTION 



The solitary kinetic Alfven wave (KAW) was studied in electron-ion (e-i) 
plasmas by Hasegawa and Mima [I] long ago. These waves can propagate in 
moderate-/? plasmas for ^ < /3(= ^0-) < 1. The phase velocity of the KAW 
is less than the thermal speed of the electrons and hot electrons are assumed 
to be inertialess and follow the Boltzmann distribution. Later on the nonlinear 
coupling of KAW and ion acoustic wave (IAW) in e-i plasma was studied by Yu 
and Shukla 

In low-/? (< ^) plasma, the electron inertia cannot be ignored and the 
wave propagation is slowed down. The phase velocity of Inertial or slow shear 
Alfven wave (SSAW) is larger than the thermal speed of the electrons. The 
nonlinear (SSAW) was studied in e-i plasma with (3 < ^ about two decades 
ago using Sagdeev potential approach. The conditions for the existence of 
localized solitary solutions were discussed. The linear dispersion relation of 
SSAW in e-i plasmas is u> — j where the external magnetic field is 

Bo = Bqz, Alfven speed is vm = ( 4An" om - )^ ( nere n io and m ; arc the ion 
unperturbed density and mass, respectively) and A e = ( here c is the speed 

of light and ui pe — ^ 4A } ^ oe J is the electron plasma frequency) is the electron 

collisionless skin depth in e-i plasmas. Furthermore k z and kj_ are the parallel 
and perpendicular wave vectors, respectively, with respect to B . 

However, the KAW has faster phase velocity than SSAW and its linear dis- 
persion relation in e-i plasmas is lo 2 — v 2 Ai k 2 (l + p 2 s k\) , where p s = ff- is the 

ion Larmor radius at electron temperature and c s — (^f - ) is the ion sound 
speed ) . The ion inertia plays an important role in shear Alfven wave dynamics 
due to ion polarization drift under low frequency assumption i.e., \dt\ << 
(where = is the ion-gyrofrequency). 

During the last decade, there has been a great deal of interest in electron- 
positron-ion (e-p-i) plasmas and a number of authors 0] EJ EJ El HH 1X0] 
have studied different linear and nonlinear wave propagation phenomenon in 
such sytems. The e-p-i plasmas are supposed to exist in active galactic nuclei 
(AGN), pulsar magnetospheres and early universe etc.pj El HI El HO E3- The 
positrons are introduced in e-i plasmas for the purpose of diagnostics and to 
model the pulsar magnetosphere in laboratory experiments [111 1121 In this 
paper, we show that the presence of positrons in e-i plasmas can change the 
nonlinear dynamics of low frequency slow shear Alfven wave significantly. Or 
we may say the other way round, that presence of ions in e-p plasmas can change 
both the spatial and temporal scales. Therefore it seems important to study the 
linear as well as nonlinear phenomenon in e-p-i plasmas to explain some aspects 
of laboratory or astrophysical e-p-i plasmas. 

Recently kinetic Alfven wave (KAW) has been studied in electron-positron- 
ion (e-p-i) plasmas by Saleem and Mahmood ^Oj • In this case, the electrons and 
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positrons have been assumed to be inertialess and they follow the Boltzmann 
distribution. The ions are considered to be inertial and cold. The phase velocity 
of the wave is less than the thermal speed of both the electrons and positrons. 
The polarization drifts of electrons and positrons can be ignored in e-p-i plasmas 
and the ion polarization drift is necessary to maintain quasineutrality conditions. 

The paper has been presented in this manner. In Sec. II nonlinear equations 
for Inertial or SSAW in e-p-i have been defined. The localized stationary solution 
in the form of energy integral equation has been obtained in Sec. III. Some of 
the possible numerical solution are presented in Sec. IV. In Sec. V discussion 
on the obtained results is presented. 

ii. SET OF EQUATIONS 

Let us consider a cold electron-positron-ion (e-p-i) plasma in the presence of 
external magnetic field B = B z along z-axis. The governing equations for the 
nonlinear Alfven wave dynamics in x-z plane in a low (3 (i.e., (3 « ^) e-p-i 

plasma with v t j < j- (where = ^j- and j = e,p) are as follows: 

The continuity equations for electrons and positrons can be written as, 

d t 7ij +d z (njV jz ) = (1) 

whereas the equations of motion for electrons and positrons along z-axis are, 

dtVjz + v jz d z v jz = —E z (2) 
rrij 

where j = e,p and qj = — e, +e for electrons and positrons, respectively. 
The ion equation of motion in the limit \d t \ « fij can be expressed as, 

v tx = -^-^-d t E x (3) 

We are ignoring the ion parallel motion along the magnetic field, so the ion 
continuity equation gives, 

d t n i + -^—d x (n i d t E x ) = (4) 
The Faradays law VxE = — ^<9 t B can be written as, 

d z E x -d x E z = --d t B v (5) 



and Ampere's law yields, 



d x B y = — (n p v pz - n e v ez ) (6) 



where the displacement current has been ignored. 
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The quasi-neutrality condition implies, 



rii ~ n e - n p (7) 

where n a (where a — e,p,i) are the densities, v az and v ax are the parallel and 
perpendicular velocities of a-species with respect to the external magnetic field, 
respectively. 

Using two potential approach, one can express the parallel and perpendicular 
electric fields as E x — —d x cf> and E z = —d z ip, where '0' and 'ip' are electrostatic 
and electromagnetic potentials, respectively. 

The linear dispersion relation of inertial Alfven wave in low-/? e-p-i plasmas 
turns out to be, 

2 _ (1+g) VAi k z fas 

-[d+p)+m (8) 

where p — (here n p o and n e Q are the unperturbed densities of positrons and 
electrons, respectively). 

The limiting case of two component e-i plasmas can be obtained by putting 
n p — (or p — 0) in above equation to obtain, 

v 2 h 2 
j = v Ai k 

which is the same linear dispersion relation of Inertial or slow shear Alfven wave 
(SSAW) in e-i plasmas which has already been studied in Ref. 2 . It can be 
seen from Eq.JHJ that the wave dispersion due to electron inertial length in the 
presence of positrons is modified because of the factor p. So the positron density 
can have significant effect on the wave dynamics in nonlinear regime as well. 

It may be mentioned here that Eq. JSJ cannot reduce to the linear dispersion 
relation of shear Alfven waves of e-p plasmas because the polarization drifts 
of electrons and positrons have been ignored in the limit \d t \ << (where 
Qj = is the gyrofrequency of j th species and j = e,p). 

in. NONLINEAR SOLUTION 

We are interested in the stationary localized planar solution of the non- 
linear set of equations. So we transform the set of Eqs. (l)-(8) in moving frame 
£ defined as , 

£ = K x x + K z z - ttt (10) 

where 'f2' is the velocity of the non-linear structure in the moving frame, K x 
and K z are the direction cosines in the x and z directions, respectively and K 2 + 
K 2 = 1. 

Now the electron continuity equation can be written as, 

v ez =u(l- n" 1 ) (11) 
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the electron momentum equation can be written by using above equation, 

K z di~ df e 2 u oe [ e ' { ' 

The normalized electron density is n e = and e z = ^f- (where m e — m p — 
m). In order to obtain Eq. lllL we have used the boundary conditions i.e., as 
£ — * | ± oo | , n e — ► 1 and v ez — > 0. 

Similarly from positron continuity equation we have, 

v pz =u(l- rip 1 ) (13) 

and the positron momentum equation and the above relation give, 

±_de^ 2 jP_ j 1 2 <P_ . t 2 

k z at oe p 2 &e [ v ' 1 j 

where u = -B- and normalized positron density is n„ = We have again 

used the boundary conditions i.e., as £ — ► | ± oo|, n p — ► 1 and v pz — > to obtain 
Eq.O- 

The ion continuity equation yields, 

de 

~ K *-j£ = (1 - nr 1 ) (15) 

where normalized ion density m = , e x = ^f- and fi e = (electron 
gyrofrequency) . 

Using Eas. i|TT|) and (|T3|) in Eq.©, one can write Ampere's law in £ co- 
ordinate as, 

K dR± = 4ne^o u ^ (j _ ^ ^ _ ^ (l _ (lg) 

Now after transforming Eq. Q in £ co-ordinate and using the above relation we 
have, 

K X K^ - KlK z ?ll = tt 2 A e - 2 [p (1 - n-/) n p - n e (l - n" 1 )] (17) 

Equating left hand side (L.H.S) of Eas. i|12|) and 114f) and then integrating the 
resulting equations twice w.r.t. £ we obtain, 

{v ez - u) 2 + (v pz - uf = 2u 2 (18) 

Using Eqs. i|ll|) and (|13|l in the above relations, we obtain, 

S 2 = 2 - n- 2 (19) 

The quasi-neutrality yields, 

1 



(1-P) 



(n e - pn p ) (20) 
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Note that the above equations hold for ^ p < 1 in three component e-p-i 
plasmas. 

Eq. (PI and Eq. © yields, 



(i-p) 



(21) 



Then Eq.jnjJ) along with Eqs.JQ) an d El gives, 
-JsT 2 [fi e fii (1 - n^ 1 )] + m 2 ^K 



— o-n„ +— — ™ 1 - n„ ) 



f> 2 A L 7 2 [p(n p -l)-(n e -l)] 



(22) 



Differentiating w.r.t £ twice the above equation and then after simplification we 
obtain, 



Kl\ 2 



if- 



x e 



a? 



3n„ 



4 /" 9nf 
3£ 



,-3' 



(i-p) ^ , _n a 2 

"2 (1 - n i ) = —2 \P n V - n e ) 



M 2 a^ 2 



where M — — is defined as Mach number. 

V A 

Integrating Ea. (|23j) twice w.r.t. £, we have, 



(23) 



3n„ 



4 f 



,<9 2 ri. 



(i-p) 

Af 2 



(l - n 4 x ) = (n p p - n e ) + (1 - p) 

(24) 



In order to obtain above equation, we have used the boundary conditions i.e., 

las £ — > | ± oo|. 

then multiplying Ea. (|2*4*jl by 'i?' both sides and 



n p — ► 1 7ij — * 1 and n. 
Let us define R = 



-3 3n e 



after integrating once w.r.t £ we obtain, 



(25) 



The Sagdeev potential is defined as, 



V{n e ) 



M 2 



M 2 



2 



(26) 



where £ = has been normalized. We have used the boundary conditions i.e. 
as £ — * | ± oo| then (^rj — * and n e — > 1 to obtain Ea. (|2*5|l . 



Equation is a well known equation in the form of "energy integral" of 
an oscillating particle of a unit mass, with velocity (^r) and position n e in 
a potential V (n e ). The conditions for the existence of localized solution of 
Eq.lJUJ require that i) V(l) = V(N Q ) = |„ e =i= 0, (N is the point where 
the curve crosses the n e axis and it can have values > or < 1) and it represents 
the maximum amplitude of the soliton h)§^r |n e =i< (where n e = 1 is the 
unstable point) and from second condition it is seen that solitary structures are 
formed only in the super Alfvenic region. 

The limiting case of two component e-i plasma can be obtained by putting 
p = in Ea. l|25(l and the Sagdeev potential in this case turn out to be, 



V(n) 



(27) 



which is the same as obtained by Shukla et. al., |2] and b — A^i^ has been 
defined. 

iv. NUMERICAL SOLUTIONS 

The numerical solutions of Eci. (|25[) are obtained for solitary structures in 
the absence as well as in the presence of positrons in e-i plasmas. The Sagdeev 
potential in the presence of positron (as a third species) in two component e-i 
plasmas becomes complicated and the third integral term does not have simple 
analytical solution. However, the numerical solutions exist and the plots of 
the Sagdeev potential V vs normalized electron densities 'n e ' corresponding 
to p = 0.2 (dashed curve) i.e., for e-p-i plasmas as well as corresponding to 
p = (solid curve) i.e., e-i plasmas for M — 1.2, K x = 0.1 have been shown 
in Fig.l. The corresponding electron density dips have been shown in Fig. 2 
for the same parameters as given in Fig.l. It can be seen from the figures 
that the amplitude and the width of solitary structures reduces with increase in 
percentage of positrons in e-i plasmas. The normalized density profiles of three 
species electron-positron-ion plasmas are shown in Fig. 3. 

The effects on the direction of propagation on the solitons for both the cases 
e-i and e-p-i plasmas are shown in figures 4 and 5, respectively. The plots for 
different directions of propagation i.e., for K x = 0.1 (solid curve) and K x = 0.3 
(dotted curve) for two component e-i ( p = 0) plasmas and for K x = 0.1 (dotted 
curve) and K x — 0.3 (solid curve) in three component e-p-i (p — 0.2) plasmas 
for the same value of M = 1.2 have been shown in figures 3 and 4. It can be seen 
from the figures that propagation direction effects on the width of the solitary 
structures. The width of the solitary structure increases with the increase in 
the obliqueness of the wave. 

The dependence of Mach number on soliton corresponding to M = 1.6 (solid 
curve) and M = 1.2 (dotted curve) for same p = 0.2, K x — 0.1 has been shown 
in Fig. 6. It can be seen from the figure that the Mach number has significant 
effect on the amplitude as well as on the width of the solitary structure. The 
wave amplitude increases and the width decreases with the increase in Mach 
number. 
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v. DISCUSSION ON RESULTS 



We have studied the solitary pulse formation of inertial or slow shear Alfven 
waves in electron-positron- ion plasmas. The conditions on the Mach number for 
the formation of such nonlinear structures in e-i plasmas were presented long 
ago • It was predicted that super Alfvenic density depletion regions can be 
formed in e-i plasmas corresponding to an arbitrary amplitude perturbation. 
However the density profiles were not plotted and the dependence of the non- 
linear structures on the propagation direction was not investigated. We have 
noticed that the pulse width increases with the increase in obliqueness (with re- 
spect to the external magnetic field) of the propagation direction in e-i plasmas 
as shown in Fig. 4. 

It is pointed out that the conditions for the formation of solitary pulses in 
e-p-i plasmas due to large amplitude SSAW perturbation are similar to the e-i 
case. That is the density depletion regions can be formed in the super Alfvenic 
region i.e, for 1 < M. But the amplitude of the soliton decreases with the 
increase in the number density of positrons. This result is very common in e-p-i 
plasmas. For example the amplitudes of the solitary ion acoustic jHJ and solitary 
kinetic Alfven waves [TO] also decrease with the increase in the concentration of 
positrons in e-i plasmas. 

In principle, the Alfven waves can have a wide range of temporal and spatial 
scales in e-p-i plasmas because the frequency of the linear wave can vary from 
uj = v Al k z (with a; < flf) tow = v Ape k z (withw < fi e , p ), wheret>A ep = ( 8 ^° m )^ 
(where m e = m p = m) is the Alfven wave speed, and correspondingly the spec- 
trum of wavelengths can be broader. These variations in the temporal and 
spatial scales depend upon the concentration ratios of different species. There- 
fore it seems interesting to analyze the linear and nonlinear wave propagation 
in such plasmas. 

Here, we have studied only the Inertial or slow shear Alfven waves in e-p-i 
plasmas which propagate on ionic time scale. It may be mentioned that the 
limiting case of e-p plasmas can not be obtained from our equations because we 
have ignored the polarization drifts of electrons and positrons. Our findings can 
be useful to explain some aspects of laboratory and astrophysical space e-p-i 
plasmas situations. 
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Figures Captions 



• Fig.l:The Sagdeev potential 'V is plotted against electron density 'n e ' 
for M — 1.2, K x = 0.1, p = (solid curve) and p = 0.2 (dotted curve). 

Fig.2:Thc normalized electron density dip decreases with the increase in 
positrons in c-p-i plasmas for the same parameters a given in Fig.l. 

Fig. 3: The normalized density profiles of electron, positron and ion are 
plotted for p = 0.2 , while M and K x corresponds to the same values as 
given in Fig.l. 

Fig.4:Effcct of propagation direction on the width of the solitary pulse in 
e-p-i plasmas is shown with K x — 0.1 (solid curve) and K x = 0.3 (dotted 
curve) for M = 1.2 and p = 0. 

Fig.5:Effcct of propagation direction on the width of the solitary pulse in 
e-p-i plasmas is shown with K x = 0.1 (dotted curve) and K x = 0.3 (solid 
curve) for M = 1.2 and p = 0.2. 

Fig.6:Thc effect of Mach number on the solitons in e-p-i plasmas is shown 
plotting n e and £ forp = 0.2, K x = 0.1, M = 1.6 (solid curve) and M = 1.2 
(dotted curve) 
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